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ABSTRACT 



X): 

Ft Context. The Herschel Open Time Key Programme TNOs are Cool: A survey of the trans -Neptunian region aims to derive physical and thermal 
, properties for a set of ~ 140 Centaurs and Trans-Neptunian Objects (TNOs), including resonant, classical, detached and scattered disk objects. One 
t goal of the project is to determine albedo and size distributions for specific classes and the overall population of TNOs. 

' Aims. We present Herschel PACS photometry of 18 Plutinos and determine sizes and albedos for these objects using thermal modeling. We analyze 
our results for correlations, draw conclusions on the Plutino size distribution, and compare to earlier results. 

Methods. Flux densities are derived from PACS mini scan-maps using specialized data reduction and photometry methods. In order to improve 
MH , the quality of our results, we combine our PACS data with existing Spitzer MIPS data where possible, and refine existing absolute magnitudes for 
|J-^ ■ the targets. The physical characterization of our sample is done using a thermal model. Uncertainties of the physical parameters are derived using 
customized Monte Carlo methods. The correlation analysis is performed using a bootstrap Spearman rank analysis. 

Results. We find the sizes of our Plutinos to range from 150 to 730 km and geometric albedos to vary between 0.04 and 0.28. The average albedo 
of the sample is 0.08 ± 0.03, which is comparable to the mean albedo of Centaurs, Jupiter Family comets and other Trans-Neptunian Objects. We 
Q ■ were able to calibrate the Plutino size scale for the first time and find the cumulative Plutino size distribution to be best fit using a cumulative power 
^ ' law with ^ = 2 at sizes ranging from 120-400 km and q = 3 al larger sizes. We revise the bulk density of 1999 TC36 and find g = 0.64^55 | j g cm"'. 
C/j , On the basis of a modified Spearman rank analysis technique our Plutino sample appears to be biased with respect to object size but unbiased with 
d t respect to albedo. Furthermore, we find biases based on geometrical aspects and color in our sample. There is qualitative evidence that icy Plutinos 
have higher albedos than the average of the sample. 
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1 . Introduction residual planetesimals during the late stages of the formation of 

. the Solar System could lead to a radial migration of the former 

^] . Smce Its discovery m 1930, Pluto has been a unique object not ^ ^.^^^j^ Neptune's outward migration, a Pluto-like body 

; only for being the only rocky planet-sized object not bound to a ^^^j^ j^^^^ ^^^^ ^ ^.^^ 2:3 resonance, excited to its 

O planet outside the orbit of Mars, but also for having the most ec- j^jgj^iy -^^^^^^ eccentric orbit and transported outward. 

CN centric and inclined or bit of all the plane t s, whi ch even overlaps i . 

^ ■ the orbit of Neptune. , Cohen & HubbardI (Il96l were the first to j— i" *e meantime, the discovery of 1992 QB 1 by| Jewitt & Luu| 

: show that, despite that overlap, close approaches between the ^1223 showed that Pluto is not the only object beyond the or- 

• . planets are prevented by the 2:3 mean motion resonance: Pluto's ^it of Neptune. The following years revealed a large popula- 

><| ■ revolution period equals 3/2 of Neptune's period, ensuring that ^lon of trans-Neptunian objects (TNOs), presumably a popula- 

H : conjunctions always occur near Pluto's aphelion. This leads to a ^lon of residual planetesimals from the age of the formationof 

P. ■ high degree of stability of the orbit. The origin of Pluto's pecu- the Solar System^^as previously proposed by | Edgeworth| (Il949|) 

har orbit (e = 0.25 and / = 17°) was first explained by MalhotTi ^nd iKuieeii dmH). The population shows dynamical complex- 

([T99l . who showed that encounters of the Jovian planets with i^y and recurring orbital characteristics allowing the classifica- 

tion of TNOs in different dynamical groups (cf. Elliot et alj 

* Herschel is an ESA space observatory with science instruments |2005|: |Gladman et al.||2p08). Many of the newly discovered ob- 

provided by European-led Principal Investigator consortia and with im- iec ts show a dyna mical behaviour similai' to that of Pluto, which 

portant participation from NASA. The Photodetector Array Camera led lJewitt&Luul (0*996) to dub them Plutinos. 

and Spectrometer (PACS) is one of Herschel's instruments. PACS has pj^j^jj^^ ^^^-^^ Neptune's 2:3 resonance dGladman et al.l 

been developed by a consortium of institutes led by MPE (Germany) ^rvrvo^ j ^ ^ i 4.- rr.- i. 

and including UVIE (Austria); KU Leuven, CSL, IMEC (Belgium); 2008) and are the most numerous resonant po pulation^ This ob- 

CEA, LAM (France); MPIA (Germany); INAF-lFSl/OAA/OAP/OAT, nervation agrees with dynamical studies by | Melita & Brunini| 

LENS, SISSA (Italy); lAC (Spain). This development has been sup- (2000), which reveals that the 2:3 resonance is much more stable 

ported by the funding agencies BMVIT (Austria), ESA-PRODEX than other resonances, particularly at low inchnations. While the 

(Belgium), CEA/CNES (France), DLR (Germany), ASI/INAF (Italy), semi-major axes of Plutinos are strongly concentrated around 

and CICYT/MCYT (Spain). 39.5 AU, their eccentricities and inclinations vary significantly 
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from Pluto's (0.03 <e< 0.88 and 0.4° < / < 40.2° with a mean 
eccentricity and inclination of 0.21 and 12.1°, respectively^ 
The origin of the Plutino population is ascribed to the same res- 
onance capt ure mechanism w hich is responsible for Pluto's pe- 
culiar orbit (lMalhotralll995b . Based on numerical simulations, 



iDuncan et al. (119951) find that the 2:3 mean motion resonance 



interferes with secular resonances leading to a severe instabil- 
ity for high inclination orbits and a longterm leakage for low 
inclination orbits of Plutinos, making t hem a possible sou rce re- 
gion of Jupiter Family Comets (JFCsl. iMorbidellil (Il997h found 
that a slow chaotic diffusion of Plutinos can be provided, which 
is necessary to explain the observed con tinuous flux of JFCs . 
Their results were recently confirmed by Idi Sisto et al ] (1201 Oh . 
who found that Plutinos may also be a secondary source of the 
Centaur popul ation. A collisional analysis of the 2:3 resonance 
carried out bv Ide Elfa et al. (2008) shows that the Plutino pop- 
ulation larger than a few kilometers in diameter is not signifi- 
cantly altered by catastrophic collisions over the age of the Solar 
System. They also pointed out the importance of specifying the 
number of Pluto-sized objects among the Plutino population, 
since the escape frequency of Plutinos strongly depends on this 
number. 

Our sample of 18 Pl utinos is selecte d based on the object 
classification scheme of iGladman et alJ ({2008) and represents 
~ 7.5% of the known Plutino population as of 27 October, 20 11 
(data from MPC). 14 of our 18 sample targets have been classi- 
fied in lGladman et al.l (2008). For the remaining ones dynamical 
simulations performed by Ch. Ejeta show that their Plutino-type 
orbits are stable for at least 10^ years, which suggests, but does 
not prove, orbital stability over the lifetime of the Solar System. 
We give the numbers and preliminary designations of the sample 
objects in Table [T] In order to avoid confusion, we will refer to 



' Statistical data are based on Minor Planet Center data as of 18 Nov. 
201 1 ( |http://minorpl anetcenter.net/iau/MPCORB.html I. 
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Fig. 1. Our Plutino sample (filled circles) compared to the known 
Plutino population (open circles, data from MPC, as of 18 Nov. 
201 1) in e-i phase space. The crosses with solid and dashed error 
bars refer to the mean and the standard deviation of our sam- 
ple and the known Plutino population, respectively. This plot 
demonstrates that our sample is heterogenous and reflects the 
dynamical spread of the complete Plutino population. (The plot 
omits 3 known Plutinos, due to their high eccentricities, which 
are not part of our sample; however, these objects are included 
in the determination of the mean.). 



Table 1. Plutino target sample. We list numbers and preUminary 
designations or object names, where allocated. In order to avoid 
confusion, we will refer to specific objects in this work only by 
their preliminary designation or name. 



(15875) 


1996 TP66 


(144897) 


2004 UXIO 


(38628) 


2000 EB 173 Huya 


(175113) 


2004PF115 


(47171) 


1999 TC36 


(208996) 


2003 AZ84 


(47932) 


2000 GN171 




2001 KD77 


(55638) 


2002 VE95 




2001 QF298 


(84719) 


2002 VR128 




2002 VU130 


(84922) 


2003 VS2 




2002 XV93 


(120216) 


2004 EW95 




2003 UT292 


(134340) 


Pluto 




2006 HJ123 



specific objects in this work only by their preliminary designa- 
tion or name, where available. In Figure [1] we plot our dataset 
in comparison to the sample of all known Plutinos in e-i phase 
space. We neglect a here, since it is very well confined around 
39.5 AU. The plot shows the spread of both samples in e-i space, 
whereas the means of the two samples show a good agreement. 
Therefore, we assume that our sample well represents the dy- 
namical variety of the whole Plutino population. 

Our Plutino sample is selected solely on the basis of Herschel 
observability in the sense of observation geometry and detector 
sensitivity. Herschel observations presented here are conducted 
over a time span of more than 7 months. In this time span, 
Herschel is able to observe at all ecliptical longitudes, despite its 
restrictions on the solar aspect angle, which must not be lower 
than 60.8° and not be larger than 1 19.2°. Hence, no bias is intro- 
duced based on preferred ecliptical longitudes. Since the target 
objects of the 'TNOs are Cool' project are selected on the basis 
of optical discoveries, the most important bias in our sample is 
presumably the optical detection bias, favoring large objects and 
objects with high albedos. The detectional bias impact on our 
results will be discussed in Sectionl4~i~ 



Our target sample includes 3 known multiple systems, which 
are Pluto, 1999 TC36 and 2003 AZ84 and flie expected binary 
2002 GN17 1 . We emphasize, that the spatial resolution of PACS 
is not sufficient to separate any multiple system. Flux densi- 
ties measured for one object therefore consist of contributions 
from fluxes of all system components. Diameter and albedo de- 
termined with our models refer to respective parameters of an 
object with the area equivalent diameter of the whole system. 
We discuss these objects in detail in Section l4!2] 

The 'TNOs are Cool: A Survey of the Trans-Neptunian 
Region' project is a Herschel Open Time Key Programme 
awarded some 370 h of Herschel observing time for in- 
vestiga ting about 140 tra ns-Neptunian objects with known 
orbits (Muller et al. '2009"). The observations include PACS 
(Poglitsch et al. 2010) and SPIRE (Griffin et al. 2010) point- 
source photometry. The goal is to characterize the diameter and 
albedo for the individual objects and the full sample using ra- 
diometric techniques, in order to probe formation and evolution 
processes in the Solar System. Using thermal modeling we deter- 
mine sizes and albedos of the 18 Plutino members. The 'TNOs 
are Cool' sample includes data of 7 more Plutinos, which have 
not been processed by the time of writing this. 
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Table 2. PACS observations summary. Column headings are: object name; Herschel ObsID of the first AOR of a sequence of 
four AORs of two consecutive observations (first visit/follow-on observation), respectively; total duration of all AORs in minutes; 
observation mid-time of all observations (MM-DD UT) in 2010; r. heliocentric distance; A: Herschel-target distance; a: solar phase 
angle; color-corrected flux density values at PACS photometer reference wavelengths 70, 100 and 160 fim, uncertainties include 
photometric and calibration uncertainties. Upper limits are 1 cr. 



Object 


ObsIDs 


Dur. 


Mid-time 


r 


A 


a 




PACS Flux Densities (mJy) 








(min) 


(UT) 


(AU) 


(AU) 


o 


70 iim 


lOOyt/m 


160 


yum 


1996 TP66 


1342202289/.. .2310 


113.3 


08-08 03:28 


27.3175 


27.6467 


2.0 


<0.8 


< 1, 


.1 


< 1.5 


1999 TC36 


1342199491/.. .9630 


75.7 


07-01 10:10 


30.6722 


30.8956 


1.9 


27.2 4 


: 1.4 


22.3 4 


: 1.9 


11.0: 


t 1.6 


2000 GN171 


1342202906/... 2971 


150.9 


08-12 18:19 


28.2876 


28.4320 


2.0 


3.2 ± 


0.7 


5.8 ± 


1.1 


3.2 4 


: 1.3 


2001 KD77 


1342205966/.. .6009 


150.9 


10-07 05:54 


35.7854 


36.1111 


1.5 


5.4 ± 


0.6 


1.0 ± 


1.1 


4.1 4 


: 1.8 


2001 QF298 


1342197661/.. .7681 


113.3 


06-03 06:43 


43.1037 


43.3727 


1.3 


7.2 ± 


0.8 


6.5 ± 


1.7 


5.0 4 


: 1.3 


2002 VE95 


1342202901/.. .2953 


113.3 


08-12 15:17 


28.5372 


28.8990 


1.9 


10.6 4 


:0.8 


8.6 ± 


1.1 


6.8 4 


: 1.6 


2002 VR128 


1342 190929/... 0990 


109.3 


02-22 11:57 


37.4636 


37.7851 


1.4 


15.8 4 


:0.9 


13.1 4 


: 1.2 


8.8 4 


: 1.3 


2002 VU 130 


1342 192762/... 2783 


112.8 


03-26 05:04 


41.6877 


42.1392 


1.2 


3.2 ± 


0.8 


2.4 ± 


1.0 


2.1 4 


: 1.3 


2002 XV93 


1342193126/.. .3175 


112.8 


03-31 22:35 


39.7152 


40.0645 


1.4 


17.3 4 


: 1.1 


17.4 4 


: 1.2 


10.8: 


±2.1 


2003 AZ84" 


1342187054 


cf. (Miiller et al. 2010) 


45.376 


44.889 


1.1 


27.0 4 


:2.8 






19.7: 


±5.2 


2003 AZ84'' 


cf . Footnote b 


484.5 


09-27 13:42 


45.3025 


45.6666 


1.2 






25.7 4 


:0.3 


14.6: 


±0.8 


2003 UT292 


1342190949/... 1025 


145.6 


02-22 22:33 


29.4217 


29.4484 


1.9 


6.3 + 


0.8 


4.7 ± 


1.4 


3.6 4 


: 1.2 


2003 VS2 


1342191937/... 1977 


75.7 


03-10 07:22 


36.4694 


36.7093 


1.5 


17.8 4 


: 1.1 


16.5 4 


: 1.5 


10.2: 


±3.0 


2003 VS2'' 


cf. Footnote c 


508.0 


08-11 18:23 


36.4760 


36.8017 


1.5 


14.4 4 


:0.3 






14.0: 


±0.6 


2004 EW95 


1342199483/.. .9712 


113.3 


07-01 22:42 


27.4708 


27.1723 


2.1 


19.5 4 


:0.9 


18.7 4 


: 1.2 


9.6 4 


: 1.8 


2004PF115 


1342208462/.. .8841 


113.3 


11-11 10:53 


41.4271 


41.2712 


1.4 


10.7 4 


:0.9 


10.6 4 


: 1.0 


8.8 4 


:2.1 


2004 UX 10 


1342199495/...9626 


75.7 


07-01 11:12 


38.9500 


39.3307 


1.4 


8.7 ± 


1.1 


10.9 4 


: 1.6 


5.2 4 


: 1.8 


2006 HJ123 


1342204150/.. .4200 


113.3 


09-09 05:52 


36.5383 


36.9867 


1.4 


3.0 ± 


1.2 


3.5 ± 


1.6 


3.2 4 


:2.1 


Huya 


1342202873/.. .2914 


75.7 


08-12 03:37 


28.6648 


28.7665 


2.0 


41.4 4 


: 1.6 


37.6 4 


: 1.8 


22.5 : 


±2.2 


Pluto/Charon 


1342191953/...1988 


75.7 


03-10 12:46 


31.7985 


32.0470 


1.7 


283.9 : 


t8.6 


354.8 4 


: 11.2 


289.2 : 


± 17.2 



Notes. Geometric data were extracted from NASA Horizons (http://ssdjpl.nasa.gov/horizons.cgi I for the indicated mid-time of the observations; 

chop/nod observations (.Miiller et al...201O) : averaged lightcurve observations using ObsIDs 1342205152/.. .5222-5225 sampling the 100 and 
160 yura bands; averaged lightcurve observations using ObsIDs 1342202371/.. .2574-2577 sampling the 70 and 160 jim bands. 



2. Observations and Data Reduction 

2.1. Herschel Observations 

Photometric measurements using the Photometer A rray Camera 
and Spectrometer (PACS, Poglitsc h et all l2010l) onboard the 
Herschel Space Observatory ( Pilbratt et al.l 1201 Ol) have been 
taken in mini scan map mode covering homogeneously a field 
of roughly 1 ' in diameter This mode turned out to be best suited 
for our needs and offers m ore sensitivity than other observation 
modes (iMiiller et al.ll2010t) . 

According to the 'TNOs are Cool!' Open Time Key 
Programme observation strategy a target is observed at two 
epochs, separated by a time interval that corresponds to a move- 
ment of 30-50" of the target, allowing for an optimal back- 
ground subtraction, to eliminate confusion noise and back- 
ground sources. At each epoch the target is observed in the 'blue' 
(nominal wavelength of 70 i-im) and 'green' (100 fim) band 
twice, using two different scan position angles. 'Red' (160 fim) 
band data are taken in parallel when sampling one of the other 
bands. This forms a series of 4 measurements in the blue and 
green bands, and a series of 8 measurements in the red band 
for a specific target. The maps are taken in the medium scan 
speed (20"/sec) mode, using a scan-leg length of 3'and 2-4 rep- 
etitions. In the case of 2002 VR128 a scan-leg length of 2.5' was 
used. More details on the observation planning can be found in 
[Vilenius et al. (submitted). 

Observational circumstances and PACS flux densities are 
summarized in Table |2l Additional information on the targets 
is given in Table[3] 



2.2. Herschel PACS Data Processing 

The raw PACS measurements were used to produce individ- 
ual scan maps (level-2) using an optimized version of th e 
standard PACS mini scan map pipeline in HIPE (Ott 20101) . 
The individual scan maps of the same epoch and band were 
mosaicked using the MosaicTask(). In the production of the 
final maps background-matching and source stacking tech- 
niques are applied to correct for the possible relative astro- 
metric uncertainties between the two visits. We created two 
sorts of background-ehminated produ cts from the 'per-visit' 
mosaics: 'background-subtracted' (s ee IStansberrv et al.l l2008t 
Santos-S anz et al.l [accepted b v A&A*) and 'double-differential' 
mapfl The latter are generated by subtracting the maps of the 
two visits, yielding a positive and a negative beam of the ob- 
ject on the differential map with all background structures elim- 
inated. A duplicate of this map is then shifted in such a way as 
to match the positive beam of the original map with the nega- 
tive one of the duplicate map. In the last step the original and 
shifted mosaics are subtracted again and averaged, resulting in 
a double-differential mosaic with a positive beam showing the 
full target flux, and two negative beams on either side showing 
half the target flux. The final photometry is performed on the 
central, positive beam. The advantage of this technique is the 
nearly complete elimination of the sky background th at makes it 
favorable for faint targets. Full details can be found in lKiss et al.l 
(in preparation!). 



- Background-subtracted mosaics were used in the cases of 
1996 TP66, 2001 KD77, 2002 VE95, 2002 XV93, 2003 AZ84, 
2003 VS2, 2004 PF115 and Pluto/Charon; double-differential mosaics 
for the other targets. 
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Table 3. Absolute magnitude H, lightcurve information, information on the presence of ices and color indices. The top part of 
the table shows: object name; a: phase angle range of observations utilized in the determination of H\ N: number of available 
observations; /3: linear phase coefficient; Hy '- absolute magnitude in the V band with respective uncertainty; photometry references; 
^mag- optical lightcurve amplitude (peak-to-peak); Ice?: information on the presence of ices, if available (discoveries in parentheses 
are tentative); LC & Ice Ref.: lightcurve and ice references. The lower part of the table shows different color indices for our sample 
objects, where available, including u ncertainties, and the s pectra l slope s in percent of reddening per 100 nm as determined from 
the given color information following lHainaut & Delsantil (l2002h . and references. The weighted sample mean (cf. Table [3]l of each 
color index is given for comparison, as well as some color indices of the Sun. 



UDject 


a 


AT 
iV 


O 

P 


Hy 


Photometry 


A 


Ice 


LL. OL ice 




( ) 




(mag/°) 


(mag) 


References 


(mag) 




Kei. 


1996 TP66 


0.90-1.86 


3 


0.10+0.04 


7.51+0.09 


1,2,3 


<0.04 


no 


4,36 


1999 TC36 


0.28-1.71 


47 


0.08+0.04 


5.41+0.10" 


1,8-13 


0.20+0.04 


H2O 


12,16,36,41,42,43 


/UUU (jJN 1/1 


A AO O A/1 

0.0/-2.04 


40 


0.14+0.03 


6.45+0.34 


10,13,14 


A 1 1 A Al 

0.61+0.03 


no 


\n TjC ta /1a /IT 

i /,36,39,40,42 


2001 KD77 


1.54-1.56 


6 


0.10+0.04 


6.42+0.08 


20,21 


<0.07 




22 


2001 QF298 


0.68-1.02 


3 


0.10+0.04 


5.43+0.07 


21,24,25 


<0.12 




22 


2002 VE95 


0.57-2.07 


40 


0.16+0.04 


5.70+0.06 


13,26 


0.08+0.04 


H2O 


26,27,36 


2002 VR 128 


0.31-0.58 


6 


0.10+0.04 


5.58+0.37 


MPC 


- 




- 


2002 VU 130 


1.28-1.37 


3 


0.10+0.04 


5.47+0.83 


MPC 


- 




- 


2002 XV93 


1.11 


3 


0.10+0.04 


5.42+0.46 


MPC 


- 




- 


2003 AZ84 


0.33-1.17 


4 


0.15+0.05 


3.74+0.08" 


14,24,28,29 


0.14+0.03 


H2O 


27,36,42,48 


2003 UT292 


0.35-1.75 


13 


0.10+0.04 


6.85+0.68* 


MPC 


- 




- 


2003 VS2 


0.53-0.59 


7 


0.10+0.04 


4.11+0.38 


MPC 


0.21+0.01 


H2O 


30,36,44 


2004 EW95 


0.45-0.80 


20 


0.10+0.04 


6.69+0.35 


MPC 


- 




- 


2004PF115 


0.31-0.93 


13 


0.10+0.04 


4.54 + 0.25* 


MPC 


- 




- 


2004 UXIO 


0.75-1.06 


4 


0.10 + 0.04 


4.75 + 0.16 


28,30 


0.08 + 0.01 


(H2O) 


30,48 


2006 HJ123 


0.22-1.49 


5 


0.10+0.04 


5.32+0.66* 


MPC 


- 




- 


Huya 


0.49-1.80 


13 


0.09+0.04 


5.14+0.07 


8,10,32 


<0.1 


(H2O) 


16,36,37,38,39,40 


Pluto/Charon 


0.18-1.87 


277 


0.037+0.002'' 


-0.67 + 0.34 


MPC 


0.33 


CH4,CO,N2 


33,45,46,47 


Object 


B-V 


V-R 


B-R 


R-I 


V-I 


B-I 


s 


Ref. 


1996 TP66 


1.03+0.11 


0.66+0.07 


1.68+0.12 


0.66+0.10 


1.29+0.11 


2.25+0.15 


30.4+4.7 


1,3,5,6,7,21 


1999 TC36 


1.00+0.13 


0.70+0.03 


1.74+0.05 


0.62+0.05 


1.30+0.13 


2.35+0.14 


32.1+2.3 


1,8,9,11-13,15 


2000 GN 171 


0.96+0.06 


0.60+0.04 


1.56+0.07 


0.62+0.05 


1.14+0.17 


2.16+0.09 


23.9+2.9 


13,14,18,19 


2001 KD77 


1.12+0.05 


0.62+0.07 


1.76+0.06 


0.57+0.07 


1.20+0.07 


2.31+0.09 


23.5+4.1 


20,21,23 


2001 QF298 


0.67+0.07 


0.39+0.06 


1.05+0.09 


0.57+0.19 


0.89+0.19 


1.35+0.09 


4.6+4.6 


19,21,24,25 


2002 VE95 


1.07+0.14 


0.72+0.05 


1.79+0.04 


0.76+0.12 


1.38+0.15 


2.47+0.13 


37.8+3.7 


13,19,28 


2002 VR128 


0.94+0.03 


0.60+0.02 


1.54+0.04 








22.76+2.1 


19 


2002 XV93 


0.72+0.02 


0.37+0.02 


1.09+0.03 








0.9+2.1 


19 


2003 AZ84 


0.67+0.05 


0.38+0.04 


1.05+0.06 


0.55+0.15 


0.92+0.14 


1.68+0.19 


3.6+3.5 


19,24,28,29,34 


2003 VS2 


0.93+0.02 


0.59+0.02 


1.52+0.03 








21.7+2.1 


19 


2004 EW95 


0.70+0.02 


0.38+0.02 


1.08+0.03 








1.7+2.1 


19 


2004 UXIO 


0.95+0.02 


0.58+0.05 


1.53+0.02 








20.2+4.4 


28,31 


Huya 


0.95+0.05 


0.57+0.09 


1.54+0.06 


0.61+0.05 


1.19+0.06 


2.14+0.07 


21.9+4.6 


7,8,10,18,32 


Sample av. 


0.84+0.13 


0.51+0.12 


1.43+0.25 


0.62+0.05 


1.19+0.11 


2.07+0.35 


17+12 




Sun 


0.64 


0.36 






0.69 






35 



Notes. We determine the weighted mean (x> of quantities j, using the absolute uncertainties cr, as weighting parameter. The uncertainty of 
(x), cr, is calculated as cr- = (l + 2i[(-"'^i" ~ (^))^/o"?]) / 2i(l/o"?)- Defined this way, cr is a combination of the weighted root mean square and 
the weighted standard deviation of the uncertainties. original photometric uncertainty is smaller than the half peak-to-peak lightcurve ampli- 
tude: new value is {a \, + (0.5A„np)^)'^^; *** converted R band data; '''' regression analysis leads to an unrealistic phase coefficient, adopt instead 
phase coefficient from lBelskava et aP l l2008h : phase coefficient adopted from Buie et al. (1997); References. MPC: photomet ric data were pro- 



vided by th e M inor Planet Center observations database (http ://mino rplanetcenter.net/iau/ECS/MPCOBS/MPCOBS.html); (1): Boehnha rdt et al. 



200 Ih : (2):|Davies et al. (2000); (3): Jewitt & Luu ( 1998); ('4'): ICollan der-Brown et al. ( 1999); (5): Tegler & Romanishin (1998); (6): Baruc ci et al. 



1999|i : (7): iJewitt & Luu ( i2o6li) : (8): iDoressoundiram et~l.l (l200l[) : (9): Delsanti et al. (2001); (10): jMcBride et al. (2003); (11): Tegleretal 



2003); ('12')~[Dotto et al.l (l2003h ; (13): Rabinowitz et al T(l2007h; (14) : DeMeo et al. (2009); (15): ^enecchi et al. (2009); (16): Ortiz etal 
2003 ): (17): Jsheppard & Jewitt (2002); (18): |Boehnhardt et al. ( 2001); (19): Tegler, private communication; (20): Doressoundira m et al. 
2OO2I) : (21): (Doressoundiram etal. (200^; (22)rTsheDDard & Jewitt (2003); (23): Peixinho et al. (2004); (24): Fornasier et al. (2004); (25): 



Doressoundiram et al. (2005); (26): Barucci et al. (2006); (27): Ortiz et al. (200^; (28):|Pema et al. (2010); (29):lSantos-Sanz et al. (2009); (30): 
Thirouin et al. (2010); (31): Romanishin et al. (2010); (32): Ferrin et al. (2001^; (33): iBuie et al. (1997); (34):'Rabinowitz et al. (2008); (35): 
Doressoundiram et al. (2008); (36): Barkume et al. (2008); (37): Licandro et al. (2001); (38): Brown (2000); (39): de Berah et al. (2004); (40): 
Alvarez-Candal et al.. (.2007): (41):, Merlin et aL(.2005.): (42):.Guilbert et aL(,2009,): r43'):.ProtoDaDa et aL(i2009i) : (44): .BaruccL tsubmitted) : (45): 
Owen et al.l dl993l) : (46) nDeMeo et alj ( 120101) ; (47): iMerlin et al.1 d20l'oF (48) nBmLCci et alj ( l201ll) ! 



Photometry is performed on the final background-subtracted eliminated. Flux densities are derived via aperture photome- 
and/or double-differential maps, which are both background- try using either IRAF/Daophot or IDL/Astrolib routines, both 
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producing identical results. We extract the flux at the photo- 
center position of the target and apply an aperture correction 
technique (Howell 1989) for each aperture radius based on the 
encircled energy f raction for the PACS point spread function 
dMiiller et al.ll201 la) . We construct an aperture-corrected curve 
of growth from which we derive the optimum synthetic aperture, 
which usually lies in the 'plateau' region of the curve of growth. 

Photometric uncertainties are estimated by random implan- 
tation of 200 artificial sources in the nearby sky background of 
the target. The 1 cr photometric uncertainty of the target flux 
is deri ved as the standard deviation of the se artificial source 
fluxes (" Santos-Sanz et al.l [accepted by A&Ah . This Icr limit is 
also given as an upper limit in the case of non-detections. 

The absolute calibration of our data is based on mini scan 
maps of 5 fiducial stars and 18 large main belt asteroids 
dMiiller et al.ll201 lal) . The absolute calibration uncertainties of 
the PACS bands are 3% for the 70 and 100 yum bands and 5% 
for the 160 fj.m band, respectively. This additional uncertainty 
adds quadratically to the photometric error and is included in 
the flux densities given in Table |2l The absolute accuracy of the 
photometry was checked against the predicted standard star flux 
densities. In the bright regime (above 30mJy) the relative ac- 
curacy was found to be ~5%, while in the faint regime (below 
30 mJy) the accuracy is driven by the uncertainties in the maps 
and is about or below 2 mJy in all bands, on the individual maps. 

A general description of the PA CS data reduction steps (in - 
cluding photometry) will be given in lKiss et alJ din preparation!) . 

The flux densities given in Table|2]have been color corrected 
and therefore are monochromatic. The measured flux density 
is determined by the response function of the PACS band fil- 
ters convol ved with the spectra l energy distribution (SED) of 
file source dMiiller et al.ll2011bl) . The SED of TNOs resembles 
a black body spectrum of a certain temperature; the resulting 
color-correction factors depend weakly on the emission tem- 
perature. The black body temperature is approximated by the 
mean surface temperature, which is given by 2 "^^T^^, where 
Tss is the subsolar temperature, defined in Equation |3] We de- 
termine the color correction factors based on the mean surface 
temperature in an iterative process during the modeling. The 
mean color correction factors for all objects in this sample are 
0.9810 + 0.0003, 0.988 + 0.003 and 1.015 ± 0.005 for the 70, 
100 and 160 i^m band, respectively: the size of the color cor- 
rections are smaller than the uncertainty of the individual flux 
density measurements (cf. Table |2]l, which justifies the use of 
the mean surface temperature as an approximation of the best fit 
black body temperature. 

We have obtained time-resolved lightcurve observations of 
2003 AZ84 and 2003 VS2, covering 1 10 % and 106% of the re- 
spective lightcurve using 99 and 100 observations, respectively. 
A detailed analysis o f both lightcurves will be subjec t of an 
upcoming publication dSantos-Sanz et al.lin preparationh . In this 
work, we make use of the averaged fluxes, which are more pre- 
cise than single measurements, due to the larger number of ob- 
servations and the cancelling of lightcurve effects. 

2.3. Spitzer MIPS Observations 

Whenever possible, the Herschel flux densities were com- 
bined with existing flux density measurements of the Multiband 
Imaging Photometer for Spitzer (MIPS, Rieke et al. (2004)) on- 
board the Spitzer Space Telescope (Wer ner et al.l 12004) to im- 
prove the results. Useful data were obtained in the MIPS 24 and 
70 yum bands, which have effective wavelengths of 23.68 and 
71.42 fim. The MIPS 70 fim band overlaps with the Herschel 



70 /im channel, providing a consistency check between data 
from the two observatories. A comparison of the PACS and 
MIPS 70 i-im band flux densities shows that most measurements 
are consistent within a 1 cr range. Significant deviations occur 
for 2003 AZ84 and Huya, which might have been caused by 
lightcurve effects or statistical noise. 

The reduction of MIPS observations of TNOs has been de_- 
scribed in detail in Stansberr y et al . (2008) and Brucker et ^ 



(2009), de tails on the calibration ca n be found in [Gordon et al 



(.20071) and lEngelbracht et all (l2007h . We adopt absolute calibra- 
tion uncertainties of 3% and 6% for the 24 and 70 fim observa- 
tions (50% larger than the uncertainties derived for observations 
of stellar calibrators), respectively. The larger uncertainties ac- 
count for effects from the sky-subtraction technique, the faint- 
ness of the targets, and uncertainties in the color corrections. 

The sky subtraction techniques introduced in Section l272l are 
derived from techniques originally developed for the MIPS data 
reductions. Reprocessed fluxes presented here are based on new 
reductions of the MIPS data, utilizing updated ephemeris posi- 
tions for the targets. This allows for more precise masking of 
the target when generating the image of the sky background, 
and for more precise placement of the photometric aperture. In 
most cases the new fluxes are very similar to the previously pub- 
lished values for any given target, but in a few cases signifi- 
cant improvements in the measured flux density and SNR were 

achieved. 

T he standard color correction routine (" Stansberry et al.l 
I2007h resulting in monochromatic flux densities for the 24 and 
70 fim bands requires the measurement of both bands and as- 
sumes the temperature of a black body which fits the 24:70 flux 
density ratio best. However, the two previously unpublished flux 
densities of 2002 VE95 and 2002 XV93 are for the 24 fim band 
only, which precludes the application of the standard method. In 
order to provide approximate color correction to these flux den- 
sities, we apply a method similar to the Herschel PACS color 
correction routine. 

2.4. Optical Photometry 

In order to derive diameter and albedo estimates, we combine 
thermal infrared measurements with optical data in the form of 
the absolute magnitude H, which is the object's magnitude if it 
was observed at 1 AC/ heliocentric distance and I AU distance 
from the observer at a phase angle a - Q. H, albedo pv and 
diameter d are related via 

d^26- lO^Q'^ ■ IQ-'^I^I ^f^, (1) 

with the magnitude of the Sun Vq - -26.76+0.02 (' BesseU et al] 
1998) in the Johnson-Cousins-Glass system and 6 a numerical 
constant. The relation returns d in km, if 6 equals 1 AU ex- 
pressed in km. In order to derive reliable albedo estimates, H and 
its uncertainty have to be well known (Harris & Harris 1997|). 
This is not the case for most small bodies in the Solar System, 
including Plutinos, for which reliable data are usually sparse. 
We determine H magnitudes from observed magnitudes in lit- 
erature and observational data from the Minor Planet Center 
(MPC) (cf. Table O. We deliberately do not use the H magni- 
tudes provided by the MPC, since there is no uncertainty esti- 
mate given for these values and t he reliability of these magni- 
tudes is questionable (for instance. iRomanishin & Teglerl (l2005l) 
have determined a systematic uncertainty of MPC H magni- 
tudes of 0.3 magnitudes). Our results include uncertainty as- 
sessments and are reasonably close other estimates (for instance 
iDoressoundiram et al.l (l2007h ). 
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Table 4. Spitzer observations summary. Column headings are: object name; Spitzer AORKEY identification; total duration of all 
AORs; observation mid-time of the 24 and 70 yuw measurements (20YY-MM-DD UT); r. heliocentric distance; A: Spitzer-target 
distance; a: solar phase angle; color-corrected flux density values at MIPS photometer reference wavelengths 23.68 and 7 1 .42 ^m, 
upper limits are 1 cr; references. 



Ohiprt 


AORKFY 


Dur. 


IVfiH-timp 




/\ 




MIPS Flux Densities (miy) 


Ref. 






(min) 


(UT) 


(AU) 


(AU) 


(°) 


24 




70 yum 




1996 TP66 


8805632 


384 


04-01-23 03:31 


26.4913 


26.2500 


2.1 


0.689 


± 0.038 


< 5.87 


1 




12659456 


40.8 


05-09-24 15:54 


26.6292 


26.1132 


1.9 


0.426 


± 0.029 


< 2.30 


1 


1999 TC36 


9039104 


62.4 


04-07-12 11:04 


31.0977 


30.9436 


1.9 


1.233 


± 0.022 


25.30 ±2.53 


1 


2000 GN171 


9027840 


25.1 


04-06-21 00:01 


28.5040 


28.0070 


1.8 


0.258 


±0.031 


5.60 ± 4.00 


2 


2002 VE95 


17766912/.. .7168 


18.0 


07-09-18 21:28 


28.2291 


28.2962 


2.1 


0.476 


± 0.044 




3 


2002 XV93 


17768704/... 8960 


227.4 


07-10-28 04:30 


40.0092 


39.6962 


1.4 


0.321 


±0.018 




3 


2003 AZ84 


10679040 


28.9 


06-03-30 09:22 


45.6674 


45.218 


1.1 


0.291 


± 0.023 


17.8 ±2.66 


1 


2003 VS2 


10680064 


11.7 


05-08-27 10:08 


36.4298 


36.5344 


1.6 


0.304 


±0.051 


25.7 ±7.34 


1 


Huya 


8808192 


21.3 


04-01-27 09:40 


29.3261 


29.2503 


1.9 


3.630 


± 0.052 


57.2 ±5.25 


1 




8937216 


43.1 


04-01-29 17:55 


29.3252 


29.2100 


1.9 


3.400 


± 0.050 


52.9 ± 1.86 


1 



Notes. Geometric data were extracted from NASA Horizons (http://ssd.jpl.nasa.gov/horizons.cgi) for the indicated mid-ti me of t he ob servations; 
References. Stansberrv et al. (2008.^ . uncertainties derived as flux density/SNR; revised 70 jim data from Stansberrv et al.l (l2008l) : previ- 
ously unpublished data. 



H and the geometric albedo are defined at zero phase angle 
{a). However, existing photometric data were taken at a 0, 
which is corrected for by assuming a linear approximation of 
the phase angle dependence: Hy - mv(l, 1) - a ■ yS, in which 
myil, 1) is the apparent V band magnitude normalized to unity 
heliocentric and geocentric distance, and /? is the linear phase 
coefficient. Since we want to derive the geometric albedo py, 
we use V band photometry when available, otherwise we use 
R band photometry and convert to V band magnitudes using 
an averaged {V - R) - 0.567 + 0.118 color inde)(0, and a cor- 
rection for the intrinsic solar color index of (V - R)q - 0.36 
dPoressoundiram et al.l 2008). {V - R) is based on obser v ations 
of 41 Plutinos and is adopted from iHainaut & Delsantil (l2002i 
data as of July 13,2010). 

If available, we give priority to data drawn from refereed 
publications, due to the existence of photometric uncertainties 
and the better calibration. All available apparent magnitudes for 
a target are converted to mv(l, 1) magnitudes and plotted as a 
function of the phase angle. For targets where the data show a 
clear trend of myi l, 1) vs. a, we use a weighted linear regression 
analysis to fit a line to the data, yielding the linear phase coef- 
ficient. In the case of high scatter, a cano nical 13 — 0.10 ± 0.04 
mag/° is assumed, based on data given bv Belskava et al.l(l2008l) 
(y and R band data, excluding Pluto). H then represents the av- 
erage of the values determined from each my(l, I). MPC photo- 
metric data are usually highly scattered due to the coarse pho- 
tometry, so the fixed-/? technique is always applied if only MPC 
data are available. The results of all computations are given in 
Table [3] which also gives lightcurve amplitude information and 
color indices, as far as available. Due to the usage of a two-visit 
observation strategy, the determined Herschel flux densities are 
a combination of flux densities taken at two different points in 
time. Hence, it is not trivial to correct both Herschel measure- 
ments for lightcurve effects. However, lightcurve effects are al- 
ready included in the uncertainty of the absolute magnitude H, 
assuming that the optical measurements are randomly sampled 
with respect to the individual lightcurve. In cases where the un- 
certainty of H is smaller than half the peak-to-peak lightcurve 
amplitude, we account for lightcurve effects by a quadratic ad- 



^ The applied Plutino V-R color agrees within 1 cr with our sample 
mean of 0.51±0.12 (cf. Tabled. 



dition of the latter value to the uncertainty of H. This was only 
necessary in three cases (cf. Table [3]). 



3. Thermal Modeling 

By combining thermal-infrared and optical data, the physical 
properties of an object can be estimated using a thermal model. 
The disk-integrated thermal emission at wavelength /I of a spher- 
ical model asteroid is calculated from the surface temperature 
distribution T{0, ip) as 

F(A) ^ed'-/A^ J J B(A, T(6, ip)) cos^ 6 cos(cp - a) dip dO, (2) 

where e is the emissivity, d the object's diameter, A its dis- 
tance from the observer, a is the solar phase angle and 9 and 
(p are surface latitude and longitude measured from the subso- 
lar point, respectively. B(A, T) is the Planck equation, B(A, T) = 
2hc^/A^ [expihc/iAksT)) - where T is the temperature, h 
Planck's constant, c the velocity of light and kg Boltzmann's 
constant. The disk-integrated thermal emission is then fitted to 
the thermal-infrared data by variation of the model surface tem- 
perature distribution T. The body's temperature distribution de- 
pends on parameters such as albedo, thermal inertia, surface 
roughness, observation geometry, spin axis orientation and ro- 
tation period. In the modeling process, we take into account the 
individual observation geometries of each flux density measure- 
ment. 

One of the first and most simple models i s the Standard 
Thermal Model (STM, cf. Lebofsky & Spencer (1989) and ref- 
erences therein), which assumes the model body to be of spher- 
ical shape, non-rotating and/or having zero thermal inertia, with 
a smooth surface and observed at phase angle a - Q. The sur- 
face temperature distribution is described by an instantaneous 
thermal equilibrium between emitted thermal radiation and ab- 
sorbed sunlight. 

In this work we ma ke use of the Near-Earth Asteroid 
Thermal Model (NEATM, Harris' ('1998)), which is, despite its 
name, applicable to a tmosphereless TNOs (Miilleretal. 2010j; 
IStansberrv et al.l2008l the latter use a similar approach with their 
hybrid STM). The main difference of the NEATM compared to 
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the STM is the use of a variable beaming factor 77, which adjusts 
the subsolar temperature 



Tss = [il-A)SQ/(€cnjr')] 



,2x1 1/4 



(3) 



where A is the Bond albedo, Sq the solar constant, cr the 
Stefan-Boltzmann constant and r the object's heliocentric dis- 
tance . A is related to the j^eome tric albedo, pv, via A ~ Ay - 
qpv dLebofsk y & Spencer"1989') . We u se a phase integral q - 
0.336 Pv + 0.479 (Brucker et al.ll2009l) . 77 accounts for thermal 
inertia and surface roughness in a first-order approximation and 
can serve as a free parameter being derived as part of the fit- 
ting process, which usually provides better results than assum- 
ing a fixed value of rj. Small values of rj {rj < 1) imply higher 
surface temperatures compared to that of a Lambertian surface, 
for instance due to surface roughness or porosity. High values 
of 77 (77 > 1) lead to a reduction of the model surface tempera- 
ture, mimicking the effect of thermal inertia. In contrast to the 
STM, the NEATM accounts for phase angles + (as shown 
in Equation O. In the case of TNOs, however, this aspect is 
unimportant, since these objects are always observed at very low 
phase angles, due to their large heliocentric distances. 

Whenever Spitzer data are available we combine them with 
PACS data. The combination of the data has a large impact on the 
reliability and accuracy of the model output. The peak emission 
wavelength in flux density units (Jy) of typical Plutino orbits 
ranges from ^Qfxm to 100 fjin (compare model SEDs in Figures 
[3]|3 and|5]l. Hence, PACS band wavelengths (70 /Jin, 100 fim 
and 160 pm) are located in the flat peak plateau and the shal- 
low Rayleigh- Jeans tail of the Planck function, respectively, and 
therefore constrain the model SED and 77 insufficiently. Adding 
a Spitzer 24 pm band measurement, which is located in the 
steeper Wien-slope, improves the ability to constrain the model 
SED and allows for the determination of reliable estimates of 
77. Therefore, we adopt floating-77 fits if additional Spitzer data 
are available and rely on floating-77 fits if only Herschel data 
are available. However, in some cases of poor data quality or 
whenever data from different parts of the lightcurve or from 
very different observation geometries are combined, the floating- 
77 method leads to 77 values which are too high or too low and 
therefore unphysical (the range of physically meaningful 77 val- 
ues is 0.6 < 77 < 2-43). In such cases, we also apply the fixed-77 
method, w hich makes use of a cano nical value of 77 = 1 .20 ± 0.35 
derived bv lStansberrv et afl (l2008h from Spitzer observations of 
a sample of Centaurs and TNOs. The validity of this approach 
is discussed in Section 14.11 We treat non-detections with a 1 o" 
uncertainty as + cr detections in our models. 

Throughout this work we assume the surface emissiv- 
ity € to be constant at a value of 0.9 for all wavelengths, 
which is based on laboratory measurements of silicate powder 
dHovis & Callahaniil966.) and a commonly adopted approxima- 
tion for small bodies in the Solar System. 

subsectionUncertainty Assessment 

In order to derive uncertainty estimates for diameter and 
albedo, we make use of a Monte Carlo (MC) simulation in which 



The physical range of 77 values were probed by using NEATM to 
model the two extreme cases of (1) a fast rotating object of high ther- 
mal inertia with low surface roughness and (2) a slow rotating object of 
low thermal inertia with high surface roughness, which yield the upper 
and lower 77 limit, respectively. Fl uxes for both cases w ere determined 
usmg a full thermophysical model jLagerr os 1996, 1997|, ll998l:lMuellej 
|200"^') . in which surface roughness is implemented via spherical cra- 
tering. The degree of surface roughness is determined via the opening 
angle (0°-180°) and surface density (0.0-1.0) of the craters. 



0) 

a 



10 



1 - 



0.1 



T r~ 

a) Original Flux Uncertainties 




Uncertainties: 

Ad =-6.8%/ +7.7% 
Apv = -16.3%/ +19.1% 
An =-12.6%/ +0.14% 



24 



70 

Wavelength (|im) 



100 



160 



0) 

a 



T 1 1 

b) Rescaled Flux Uncertainties (S = 3.33) 




70 

Wavelength (|j.m) 

Fig. 2. Illustration of the rescaled flux density uncertainties tech- 
nique using synthetic data, a) Filled circles represent the flux 
densities at the different wavelengths with original flux density 
uncertainties. The dark line represents the best fit line (x^^d ~ 
11. 09); the grey lines are fit lines of 100 out of 1000 synthetical 
bodies generated for the MC method. The bad quality of the fit is 
obvious, since the 70 and 100 pm datapoints are hardly fit by any 
of the single fit lines. Therefore, the reliability of the uncertain- 
ties of the output parameters is questionable, b) By rescaling 
the error bars by a factor of 3.33 = Vl 1 .09, the broadness of 
the distribution of single fit lines is increased, which covers now 
all datapoints. As a result, the output parameter fractional uncer- 
tainties are increased and therefore more realistic. 



a sample of 1000 randomized synthetic bodies is created by vari- 
ation of the observed flux density, the 77 value, and the absolute 
magnitude H. The randomized parameters follow a normal dis- 
tribution centered on the nominal value, within the limits of the 
respective uncertainties. The uncertainties in flux density and H 
are taken from Table |2] and Table [3] respectively. In the case 
of floating-77 fits, 77 follows as a result of the modeling process, 
whereas i n the case of fi xed-77 fits, the value of 77 + 0.35, deter- 
mined bv lStansberrv et al. (2008) is applied. 

Uncertainties of a specific parameter are derived from the en- 
semble of modeled synthetic bodies as the upper and lower val- 
ues that include 68.2% of the ensemble (centered on the median), 
respectively. The uncertainties are usually aligned asymmetri- 
cally around the median, since most diameter and albedo distri- 
butions do not stri ctly follow a n ormal distribution. This method 
was introduced bv lMueUer et al.l (1201 li) . Finally, we combine the 
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uncertainties determined through the MC method with the re- 
spective parameter values provided by the best model fit. 

In order to improve uncertainty estimates, we make 
use of the resca l ed uncertainties approach, introduced by 
ISantos-Sanz et alj ([accepted by A&Xb . The motivation for this 
technique is the fact, that uncertainty estimates resulting from 
the MC method might be unrealistically small in cases in which 
the model fits the observational data only poorly. As a measure 
for the fit quality we use the reduced which is a result of the 
fitting process. It is defined as 



2 

^red 



1 



n.d.f. 



^obs,i 



(4) 



with Fohs^i and F„,ociei,i being the observed and modeled flux den- 
sities of the different observations /, respectively, and cTobsj be- 
ing the respective observational flux density uncertainty, x^^^^ is 
normalized by the number of degrees of freedom (n.d.f.), which 
equals the number of datapoints -2 in the case of a floating-;; fit 
and the number of datapoints -1 in the case of a fixed-?/ method, 
respectively. A poor fit results in » 1 . 

The idea is to increase the flux density uncertainties uni- 
formly by a factor S, which leads to an increase in variation 
of the flux densities of the synthetic bodies. The value of S is 
determined in such a way as to provide by definition the best 
possible model fit, i.e. x^^^ - 1, which is achieved by choosing 



S - Jy I ■ , with Y . . being the minimum v , resulting 

\^ i-ed,mtn^ ^ red,mm ^ red 

from fitting the unaltered observational data. Hence, the rescaled 
flux density uncertainties are &obs = S (Tobs and the resulting 
x]ed,mm^ which is obtained by replacing cTobs with &obs in equa- 
tion |4] is by definition unity. Diameter uncertainties are dom- 
inated by flux uncertainties and the quality of the fit, whereas 
albedo uncertainties are dominated by the uncertainty of H. The 
rescaled uncertainties are used only in the MC method to de- 
termine the uncertainties of the physical parameters. In order to 
obtain the best fit, the original uncertainties are applied. 

The method is illustrated in Figure |2] using synthetic data. 
Figure |2] a) shows the original flux densities with error bars of 
the size of the original uncertainties. The dark line represents 
the best NEATM fit (floating-?; method) of the dataset, whereas 
each grey line symbolizes the best fit SED of one of 100 of the 
1000 randomized synthetic bodies. It is obvious that the distri- 
bution of single fit curves misses the 100 /im PACS band com- 
pletely and hardly fits the 70 iim datapoint. Hence, the fit quality 
is bad, which results in an underestimation of the model output 
uncertainties. By rescaling the uncertainties of the flux densities, 
the broadness of the distribution of single fit lines increases, as 
depicted in Figure |2]b). As a consequence of the rescaling, all 
datapoints are covered by the set of lines and the derived uncer- 
tainties are significantly increased, leading to much more realis- 
tic uncertainty estimates. 

4. Results and Discussion 

An overview of all modeling results is given in Table |5] Where 
possible, floating and fixed-/; modeling approaches were applied 
based on Herschel-only and combined Herschel and Spitzer 
datasets, in order to test the consistency of the results. The ob- 
jects' SEDs determined by the different models are plotted in 
Figures [m and 121 

In the following sections we test the validity of the 77 = 1 .2 
assumption for the Plutino sample, and discuss the modeling re- 
sults of the individual objects. After that, the results of the sam- 



Table 5. Modeling results sorted by modeling technique. 'Data' 
denotes the dataset on which the modeling results are based: 'H' 
equals Herschel data, 'HS' equals Herschel and Spitzer data, re- 
spectively; diameter, geometric albedo and rj are given including 
respective upper and lower uncertainties. We show all succes- 
ful model fit results in order to allow for a comparison of the 
reliability of the techniques and the data sets. 



Object 


Data 


d [km] 


Pv 


n 


floating-;; fits: 










1999 TC36 


H8 


393.1+;" 

— iO.JS 


0.079+»-!!!3 

— U.UI 1 


l-10to8 

— U.Ufi 


2000 GN171 


HS 




215+''-''''' 

-0 070 


1 11+0-2'' 
-0.21 


2002 VE95 


HS 


249.8!|^;^ 


0.149!»«« 


1 40+°'2 


2002 XV93 


H 


451 3+"'-' 


o.o60!S:S!g 


72+°-^'' 

"• ' "^-0.24 


2002 XV93 


HS 


549.2!i^;^ 


040+""'" 

"■"^"-0.015 


1 94+0.06 
^- -0.06 


2003 AZ84° 


HS 


727.0!^^;^ 


107+°''^^ 

"■'^"'-0.016 


1 05+°''' 
^■"-'-0.15 


2003 VS2° 


HS 


523.0!34i 


147+0.063 
"■'-^'-0.043 


1 «7+0.30 
-0.23 


2004 EW95 


H 


265.9!^^? 


053+""'^ 


0.93!°:^^ 


2004 PF 115 


H 


406.3!^?;^ 


113+""**2 

"■^-'-^-0.042 


0-84!°:5i 


2004 UX 10 


H 


361.0!^f/ 


172+°''*' 


96+° '*' 

"■=^"-0.56 


2006 HJ123 


H 


283.1!i«;^ 


136+°-'°'* 


2 48+'-''" 

^•^°-1.91 


Huya 


H 


395.7+3|^ 




75+°-22 


Huya 


HS 




o.o8it«!l 


89+°"'' 

"■°^-0.06 


fixed-/; fits: 










1996 TP66 


HS 


154.0!^8-8 


0.074+00" 

— U.U J 1 




1999 TC36 


H 


407.5i^-' 


0-073+!;!!!' 

—U.U 1 4 




1999 TC36 


HS 


428.3!?|2^ 


o.o66!;!:»g 




2000 GN171 


H 


150.0+;«? 


0.207+"-'°^ 

— U.U/U 




2000 GN171 


HS 


155.11 


0.193+!!-l^' 

—\)XiIj 




2001 KD77 


H 


232.3+t[!i 


0.089+!!!!:;:! 

— U.U2 / 




2001 QF298 


H 


408.2!^2 


0.071!":™ 




2002 VE95 


H 


237.0!^?;2 


165+°"" 




2002 VE95 


HS 


226.3!^^6 


0. 181^QQ^2 




2002 VR128 


H 


448.5!^!^ 


0-052!»-»?^ 




2002 VU130 


H 


252.9!^3-6 


0.179!»-202 




2002 XV93 


H 


544. r^?-^ 


0.041+!!"2^ 




2002 XV93 


HS 


5343+69^6^ 


0.042+o"^8 

-0.021 




2003 AZ84° 


H 


771.0!^"^;^" 


0.095+!!"26 

— u.uiy 




2003 AZ84° 


HS 


766.8!f,^;^ 


096+°"^^ 

"■"^"-0.016 




2003 UT292 


H 


186.6!i^3 


092+""'" 

"■"'^^-0.049 




2003 VS2 


H 


483.1!^^;° 


n 179+0.092 
"■^ '^-0.056 




2003 ¥82" 


H 


465.7!^^;^ 


0.185t™ 




2003 ¥82" 


HS 


466.9!^8-6 


0.184^- 




2004 EW95 


H 


291 1+2"^ 


044+''"2' 

"■"^ -0.015 




2004 PF 115 


H 


468.2!38-6 


123+""''^ 




2004 UX 10 


H 


398. l!i^ 


0-141!°S 




2006 HJ123 


H 


216.4!^^;^ 


0.281!-- 




Huya 


H 


461.111:^ 


073+°"''' 
"■"'-'-0.011 




Huya 


HS 


561 5+**^'* 


049+°"^*' 

"■"^^-0.017 




Pluto/Charon 


H 


71 1Q Q+164.5 


730+°"'^ 

'-'"-0.153 





Notes. In the case of all fixed-// modeling approaches, /; 
was applied. Using averaged lightcurve observations. 



1.20 ±0.35 



pie as a whole are discussed and used to determine a cumulative 
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size distribution of the Plutino population. Finally, correlations 
between the different physical, photometric and dynamical pa- 
rameters and the presence of ices are examined. 

4. 1 . Test of the Fixed-r] Assumption 

We were able to derive floating-77 fits for 12 targets; we adopted 
the 7 which are based on combined Herschel and Spitzer data. 
The top section of Table |5] gives the results of all converged 
floating-77 fits. The weighted mear0 rj determined from the 7 
adopted floating-/; model fits yields (77) = which 
is consistent with the adopted value of 77 = 1.20 + 0.35 
dStansberrv et al.l2008l) applied in the fixed-77 approach. The root 
mean square fractional residuals in the results of the fixed-77 
method between using 77 - 1.20 and (77) = 1.11 is 5 % and 
10 % for diameter and albedo, respectively. However, the aver- 
age fractional uncertainties in diameter and albedo of all adopted 
results are + 10 % and +49/-31 %, respectively, which is larger 
than the residuals. Therefore, we continue to adopt 77 - 1.20 for 
the fixed-77 approach. 

4.2. Discussion of Individual Targets 

1996 TP66: Despite the availability of both Herschel and Spitzer 
data, the floating-77 approach was not applicable for 1996 TP66, 
because of the non-detections in the PACS 70, 100 and 160 fim 
bands. The best fit fixed-77 solution is largely constrained by 
the two 24 fim MIPS flux densities. Therefore, our results 
(d = 154.0!|;^ km , pv = 0.074+°°^^^ bas ed on 77 = 1.20) 
barely differ from those of IStansberrv et all ([2008), who find 
d = 160^45 km and pv = 0.074!° using 77 = 1.20t[J-^^^ solely 
based on the 24 jjin MIPS flux densities. The reason for the non- 
detections is unclear. However, the fact that no clear detection 
was possible in the wavelength regime of 70 pm and above us- 
ing two different instruments (PACS and MIPS) and at different 
epochs points to a real effect. We based the observation planning 
on Pv - 0.07, which agrees well with the best fit model. 
7999 TC36: This target was weU detected in all three PACS 
bands. In combination with MIPS flux densities, a good 
floating-77 fit was possible here, which misses the PACS 
160 pm band flux density. This might hint to a more com- 
plex surface albedo distribution or a wavelength-dependency 
of the emissivity. Our adopted diameter and albedo estimates, 
393.1;^26 8 a nd 0.079!j|;j|}|, re spectively, agree with ear- 
lier estimates by IStansberrv et al.l (l2008i) : d - 414.6+3^ 2 km 
and Pv ^ 0.072+g;g}^. 1999 TC36 is a triple compo- 
nent system ('Truiillo & Br ownll2002l:|jacobson & Margotll2007t 
[Benecchi et al. 2010), consisting of two similarly sized central 
components Ai and A2 and a more distant secondary compo- 
nent, B. The nature of this system makes it possible to determine 
the mean density of the system and, using optical flux differences 
between the single components, the sizes of the components. We 
revisit the calculations performed by Stansberrv et al. (2006) and 
iBenecchi et al.l (l2010l) . Our newly derived estimates of the sizes , 
based on the flux differences found by iBenecchi et al.l (1201 Oh . 
are: dAi = 272+|^ km, dA2 = 251![^ km and = 132+^ km. 
Using these diameters, we are able to improve the mean system 
density g - 0.64!y |j g cm""*, which is somewhat higher than 
the earlier estimate g - 0.54- l ii g cm~^ (IBenecchi et al ] l2010l) . 
Assuming material densities of 1.0 < go < 2.0 gem"-', we find 

^ The weighted mean and its uncertainties were determined using the 
equation given in Table [3] 



porosities or fractional void space of 36-68%. Our calculations 
are based on the assumptions of spherical shape of each compo- 
nent and the same albedo for each component. 

2000 GN171: Both fixed and floating-77 model fits based on the 
combined Herschel and Spitzer datasets were possible, leading 
to good fits with similar outcomes. We adopt the floating-77 result 
d = 147.1+207 ^ 0.215![;;[;^,^, and 77 = l.ll+»-^|, which 

differs significantly from an earlier estimate bv IStansberrv et afl 
(2008), who determined d = 321+5^-4km, pv = 0.057+°°25 and 
77 - 2.32+Q43. The large discrepancy stems f rom an overesti- 
mated MIPS 70 pm band flux density in the IStansberrv et al.l 
(I2OO8I) data, which has been revised for this work. By adding 
the Herschel PACS data, the SED is better constrained, leading 
to a more reliable result. 2000 GN171 is either a Roche binary 
or a Jacobi ellipsoid (Lacerda & Jewitt 2007). With a lightcurve 
amplitude of 0.61 mag ( Sheppard & Jewitt,2002 ). the size of the 
companion or the ellipticity of the body is significant. Hence, the 
real diameter of 2000 GN 171, the smallest Plutino in our sample, 
might be even smaller. 

2007 KD77: We modeled this target solely on the basis of PACS 
data, since existing MIPS flux density measurements suffer from 
a large positional uncertainty. For this faint target, the fixed-77 fit 
mainly relies on the 70 pm band measurement, which has the 
highest SNR of all bands. Figure [3] clearly shows the mismatch 
of the 100 pm band measurement, which is scarcely detectable 
on the sky subtracted maps. Due to its large uncertainty, this 
datapoint is nearly neglected by the fixed-77 fit, and therefore 
scarcely impacts the modeling results. The reason for the low 
100 pm flux density is unknown. 

2002 VE95, 2002 XV93, 2003 AZ84, 2003 VS2 and Huya: These 
objects show a good to excellent agreement of both the differ- 
ent model technique fits compared to each other, as well as the 
model SED compared with the different measured flux densities. 
For all five targets combined Herschel and Spitzer flux densi- 
ties are available and we adopt the floating-77 modeling results. 

2003 AZ84 and 2003 VS2 data include averaged lightcurve data, 
which strongly constrain the respective SEDs, due to thei r small 
uncertainties. Earlier estimates of 2003 VS2 by Stans berrv et al.l 
(I2OO8I) based on Spitzer data suggested d = 725.2+J^^ ° km and 
Pv = 0.058+° °^^ using 77 = 2.00+[5-^|, which is significantly dif- 
ferent from our results d = 523.0+i^^-^ km, pv = 0.147+°°g, 
based on 77 = 1. 57+023. The discrepancy might stem from 
the high 70 pm MIPS flux density value, which, together with 
the 24 pm MIPS data point, suggests a steeper Wien slope 
of the spectral energy distribution compared to the combined 
Herschel and Spitzer data set. 2003 AZ84 was observed before 
by Herschel during t he Science Demo nstration Phase using the 
chop/nod technique (Muller et al]|2010 ). yielding d - 896 ± 55 
and Pv = 0.065 ± 0.008, using 77 = 1.31 + 0.08, which differs 
from our results d = 727.0!^^-^ km, pv = 0.107+JJ JJ23, based 
on 77 = 1.05+Q [g. Despite the smaller diameter determined us- 
ing the averaged lightcurve data, 2003 AZ84 is still the largest 
Plutino in our sample, apart from Pluto. 2003 AZ 84 has a moon, 
which is 5.0 ± 0.3 mag fainter than the primary (iBrown & Sueii 
120071) . This large difference in magnitude suggests that the ther- 
mal flux of the companion is negligible. The modeling results 
of Huya vary significantly depending on the model technique 
and the dataset. A visual inspection of the model fits shows 
that the floating-77 approach based on the combined data set 
matches the measured SED best and results in li = 438.7+252 
and Pv ^ 0.081+°°|], b a sed o n 77 = 0.89![J|J^. Earlier esti- 
mates by IStansberrv et al.l (l2008l) based on Spitzer data suggest 
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Fig. 3. Model fits of the Plutino sample using different modeling approaches and different data samples: the solid and dashed lines 
are floating-// fits based on combined Herschel and Spitzer and Herschel-only data, respectively; the dash-dotted and dash-dot-dotted 
lines represent fixed-// fits of combined Herschel and Spitzer and Herschel-only data samples, respectively; datapoints at 23.68 fim 
and 71.42 fxm are Spitzer MIPS data (open circles), fliose at 70 fim, 100 fim and 160 fim are Herschel PACS data (fiUed circles). The 
model is fit to each flux measurement using the circumstances appropriate for the epoch of that observation. In order to simplify the 
figures, we have normalized the Spitzer flux densities to the epoch of the Herschel observations using the ratio of the best-fit model 
flux densities for the Herschel and Spitzer epoch, and then plotted the measured value times this ratio. Uncertainties are rescaled 
similarly, preserving the SNR for the Spitzer data. 



d = 532.6+251 and pv = 0.050+°°°^, based on // = 1.09:^g-g^. 2007 QF298, 2002 VR128, 2002 VU130, 2003 UT292, 
2002 XV93 has the lowest albedo in flie sample. ' 2004 EW95, 2004 PF115, 2004 UXIO and 2006 HJ123: For 
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Fig. 4. Continuation of Figure [3j 



these targets we have to rely on Herschel-only data and there- 
fore adopt the results of the fixed-;; technique. We were able 
to succesfully apply the floating-77 method to 2004 EW95, 
2004 PF115, 2004 UXIO and 2006 HJ123 as well. Modeling 
results in Table |5] and the plots in Figures [3] and |4] of the SEDs 
show a good agreement between the results of the different tech- 
niques, supporting the validity of the fixed-77 approach and the 
reliability of the Herschel flux measurements. Differences in the 
fit quality are induced by the different SNR values of the respec- 



tive measurements. Both 2002 VU130 and 2006 HJ123 suffer 
from large albedo uncertainties, which is a result of the large un- 
certainty of their absolute magnitude H. 2006 HJ123 shows the 
highest albedo of our sample, excluding Pluto. 
Pluto: Paradoxically, Pluto is a clear outlier from the Plutino 
population, being by far the largest and having the brightest sur- 
face. Besides th at, Pluto was the first TNO discov ered to have 
an atmosphere (iHubbard et al.lll988l: lBrosch|[T995h . has a pro- 
nounced optical lightcurve caused by albedo variations, and has 
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Fig. 5. Continuation of FigurelH 



a moon half its own size. An additional complication is the ef- 
fect of N2 ice on Pluto's surface: sublimation and deposition of 
N2 is accompanied by latent heat transport that results in much 
lower (higher) day side (nightside) temperatures than would be 
expected on an airless body (e.g. Spencer et aLlll997h . Despite 
these facts, we attempt to treat Pluto like any other Plutino 
and compare our modeling results to the known properties of 
Pluto and Charon, and discuss the reasons for discrepancies that 
emerge. 

The Pluto system has been subject to detailed study via, 
for instance, direct imaging, stellar occupation observations and 
thermal radiometry in the past. Mean hemispheric albedo ranges 
are 0. 49 to 0.66 for Pluto and 0.36 to 0.39 for Charon (Buie et al. 
Il997h. Charon's diameter is well determined to be 1208 km 
jGulbis etalll2006l: ISicardv et al]|2006l) . Pluto's diameter esti- 
mates suffer from uncertainties related to the pre sence of the at- 
mosp here and range from about 2290-2400 km (iTholen & Buid 
For the purpose of comparing with the r esults of our 
therm al models, we adopt a diameter of 2322 km (lYoung et al.l 
120071) . These numbers lead to an equivalent diameter and albedo 
of the combined system of 2617 km and 0.44-0.73, respectively. 

We deliberately refrain from combining our Herschel 
data set with exist ing Spitzer measurements (for instance 
iLellouch et al.l(l201 due to th e so far unexplained se cular fad- 
ing of the Pluto/Charon system (IStansberrv et al.ll2009h . Instead, 
we rely on our original Herschel measurements only, and so 
we adopt the fixed-77 results: d = 2119.9^[g2 2 km and pv - 
0.730;^Q [g3 (cf. Figure |5]l. Our thermal estimate of the effective 
system diameter deviates by 3 cr from the known effective sys- 
tem diameter given above; the albedo derived from our thermal 
modeling is only barely consistent with the maximum allow- 
able albedo given above. Our attempt to apply a simple thermal 
model to Pluto shows the susceptibility of such models to viola- 
tions of the model assumptions. We discuss effects, which may 
explain the discrepancy of our adopted values with the equiva- 
lent values: 

- By replacing the Pluto/Charon system with a single equiv- 
alent body, we had to assume a single average albedo for 
both objects, which is a coarse simplification, in order to 
apply our models. Furthermore, surface albedo variations 
of Pluto (IBuie et al.ll2010l) lead to a more complex spec- 
tral energy distribution and affect the model results signif- 
icantly. Our model results would be improved by using a 
two-terrain model, whic h was us ed to model the thermal 
emission of Makemake (iLim et a l. 2010), in which Pluto 
and Charon would each be represented by one terrain, or 
applying more sophisticated albedo surface distributions: 



ILellouch et al.l ( 1201 1 ) applied a 3 -terrain model to model the 
thermal emission of Pluto and an additional terrain to model 
Charon's contribution. 

- Former thermal infrared measurements of Pluto revealed 

its thermal lightcur ve in different wavelength ranges 

(ILellouch et al.ll2000l lIoTll) . which was not accounted for 
in our modeling. Given the time lags between the individual 
observations, this might have introduced a perceptible dis- 
tortion of the SED, which affects the modeling results. 

- Millimeter wavelength observation s of Pluto 

jCiurwell & Bulled 120051: iGurwell et alJ 1201 Oh revealed 
a lower surface brightness temperature than assumed from 
thermal equilibrium, which is due to the combination of two 
effects: (1) the mixing of SEDs at different temperatures, and 
in particular the increased contribution of the N2-dominated 
areas, whose temperature is ~ 37 K due to the dominance 
of sublimation cooling in the thermal budget of N2 ice, 
and (2) the fact that emissivities are expected to decrease 
at long wavelengths as a result of sub-surfa ce sounding in 
a prog ressively more transparent medium. ILellouch et al.l 
(I2OI I t) find in particular that the emissivity of CH4 ice, 
one of the terrains in their 3 -terrain Pluto model, shows a 
steady decrease with increasing wavelength (their Figure 8). 
Without additional assumptions, our thermal model is not 
capable of taking such effects into account. 

Due to the failed model fit and Pluto's uniqueness among 
the Plutinos, we exclude it from the following discussion of the 
properties of our Plutino sample. 

4.3. Sample Statistics 

Figure |6] shows diameter and albedo probability density func- 
tions and histograms based on the results compiled in Table 
|6] The size distribution ranges from 150 km to 730 km and 
the albedo distribution from 0.04 to 0.28 with a clear peak 
around 0.08. The low probability density at high albedos is 
caused by the large uncertainties of such albedos. Excluding 
Pluto, the weighted meaij^ of the albedo, weighted by abso- 
lute uncertainties, yields 0.08 ± 0.03, which agrees well with 
the range of typical geometric albedos of the TNO and Centaur 



* Upper and lower uncertainties were calculated corresponding to the 
calculations applied in Table [3] The calculation of the weighted mean 
albedo based on fractional uncertainties yields 0.11. However, we prefer 
to adopt the weighted mean albedo based on absolute errors, since it 
better agrees with the low albedo of the bulk of the Plutinos and is 
enforced from our results in Section l4!4l 
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Fig. 6. Diameter and albedo probability density functions and 
histograms. The probability density functions describe the prob- 
ability to find an object of a specific diameter or albedo, when 
picking a random object. The functions were derived by sum- 
ming up the probability density functions of each object, being 
represented by a non-symmetric normal distribution taking into 
account the different upper and lower uncertainties of each re- 
sult. This approach returns a more detailed picture than the his- 
togram. Pluto was excluded from the sample. 

sample of IStansberrv et al.l (l2008h . Table |7] shows a compar- 
ison of the mean albedos of our Plutino sample with differ- 
ent TNO sub-population and JFCs. It turns out that the aver- 
age albedo of the Plutinos agrees well with that of the scat- 
tered disk. Centaur, hot classical and JFC population. Through 

Table 6. Adopted modeling results. Columns as in Table |5] rj 
values of floating-/; modeling results are in bold. 



Object Data d [km] pv r] 



1996 TP66 


HS 


154.0!i:^ 


074+''-''63 


1 20+" 35 

^■^ -0.35 


1999 TC36 


HS 


393.1!^^:^ 


079+°-'"3 

"■"'^-0.011 


l_lQ+0.07 

-0.08 


2000 GN171 


HS 




"■^^-^-0.070 


111+0.24 
-0.21 


2001 KD77 


H 






1 20+" 35 
'-■^ -0.35 


2001 QF298 


H 




071+"-™ 


1 20+" 35 
^-■^"-0.35 


2002 VE95 


HS 


249.8+;;^;^ 


149+0.019 

"■^^^-0.016 


140+0.12 

-0.11 


2002 VR128 


H 


448.5!«:^ 


052+"-"-^ 

"■"-'^-0.018 


1 20+" 35 
^-■^"-0.35 


2002 VU130 


H 


252.9+33-6 


179+0.202 
"■'^ '^-0.103 


1 20+" 35 
'-■^ -0.35 


2002 XV93 


HS 


549.2!^^:^ 




1 24+0.06 

^•^^-0.06 


2003 AZW 


HS 


727.0!^^:? 


107+''"^3 

"•^"'-0.016 


l.OS+o" 

-O.IS 


2003 UT292 


H 


185.6!!^:^ 


067+''-''6'* 

"■"-"-"-0.034 


1 20+" 35 
'-■^ -0.35 


2003 VS2° 


HS 


523.0!35-i 


147+0.063 
"■'-^'-0.043 


157+0.30 

-0.23 


2004 EW95 


H 


291.1!^o-3 


"■"^^-0.015 


1 20+" 35 


2004 PF 11 5 


H 


468.2+38-6 


1 7 "3 +0.043 
U.1Z.J_QQ33 


1 20+" 35 
^-■^"-0.35 


2004 UXIO 


H 


398.i:32-« 


0.141!-- 


1 20+" 35 
^-■^"-0.35 


2006 HJ123 


H 


216.4!^^ 


0.281!«:^=^ 


1 20+" 35 
^■"^"-0.35 


Huya 


HS 




0.081!-|| 


0.89+"'"' 

-0.06 



Notes. No adopted values are given for Pluto/Charon, since our model 
approaches are not applicable to these objects (cf. Section l42t . using 
averaged lightcurve observations. 



Table 7. Average albedos of Centaurs, different TNO sub- 
populations and Jupiter Family Comets. 



Population 


Average Albedo 


N 


Reference 


Centaurs 


0.06 


21 


1 


Plutinos 


0.08 ± 0.03 


17 


2 


Cold Classicals" 


0.16 ±0.05 


6 


3 


Hot Classicals" 


0.09 ± 0.05 


12 


3 


Scattered Disk Objects* 


0.07/0.05 


8 


4 


Detached Objects* 


0.17/0.12'- 


6' 


4 


TNOs and Centaurs 


0.07 - 0.08 


47 


1 


JFCs 


< 0.075'* 


32 


5 



Notes. N denotes the number of datapoints. unweighted mean; *** un- 
weighted/weighted mean, based on relative uncertainties; excluding 
Eris; means the geometric albedo in R ba nd, which is compara - 
ble to Pv for small values. References. (1): Stansberrv et al. ('200^ 
(2): this work; (3) : [Vilenius et al. (submitted); (4): Santos-Sanz et aP 
jaccepted bv A&Al) ; (5): IFemandez et alj ( 120051) . 



dynamical studies 'Dun can et al] (Il995h . lYu & Tremaind (Il999h 
and di Sisto et al, suggested the Plutino population to be 

a source of JFCs and/or Centaurs, which is compatible with our 
results. Furthermore, most of the Plutinos appear to have darker 
surfaces than the cold classical TNOs and detached objects. 

4.4. Plutino Size Distribution 

Using the absolute magnitude estimates for all known Plutinos 
from the MPC and our measured albedos we are able to calibrate 
the Plutino size scale for the first time and to determine a cu- 
mulative size distribution of the known Plutino population. We 
use two different approaches to determine the cumulative size 
distribution N{> d). Firstly, we assume monochromatic albedo 
distributions based on the two average albedos determined in 
Section 14.31 This approach enables the direct conversion of H 
magnitudes into diameter but neglects the measured diversity 
in albedo. The second approach makes use of a Monte Carlo 
method to determine the size distribution based on the actually 
measured albedo distribution. For this purpose, each MPC H 
magnitude is assigned to a randomly generated albedo, which 
is used to determine the diameter Hence, we assume that H is 
not correlated to the albedo, which is supported by our findings 
in Section |43] The distribution of the randomly generated albe- 
dos follows the determined probability density function (Figure 
|6]l. Hence, this approach takes into account the measured albedo 
diversity. FigureQshows a comparison of the two different meth- 
ods. It is clearly visible that the averaged Monte Carlo approach 
agrees well with the monochromatic py = 0.08 model, which 
supports our assumption that this value better represents the av- 
erage Plutino albedo compared to the value of 0. 11 . 

For comparison, we plot in Figure |7] the slopes of different 
cumulative power laws N{> d) ~ d^''. We find that the dis- 
tribution of intermediate-sized Plutinos (120-450 km) is well- 
described by (7 = 2, wh ich is smaller than the presumed value 
of ^ ~ 3 dTruiillo et a l. 2001; Keny on & Bromley 2004). The 
large diameter tail (d > 450 km) seems to be better described 
by q - 3. However, this region suffers from uncertainties due 
to small number statistics. A change in the cumulative power 
law slope at the small-sized tail seems to occur at diameters 
of 80-1 20 km, which i s larger th an the p roposed ra nge of 40- 
80 km (iBernstein et al.l l2004; Kenvon et al . 2008; de Ella et al.1 
I2008h . However, the number of intrinsically fainter, and there- 
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Fig. 7. Cumulative Plutino size distributions A^(> d) based on 
MPC H magnitudes and different assumptions: the solid and 
dashed line represent monochromatic albedo distributions using 
Pv = 0.08 and pv = 0.11, respectively. The grey line shows 
the result of an averaged Monte Carlo simulations based on 100 
runs and the measured albedo distribution with corresponding 
1 cr standard deviations. It is obvious that the monochromatic 
Pv - 0.08 size distribution agrees better with the result of the 
Monte Carlo simulation. Furthermore, we show two cumulative 
power laws A^(> d) ~ d^'' with q = 2 and q = 3 and arbitrary in- 
terceptions for comparison. Aq = 2 cumulative power law better 
fits the measured size distribution in the size range 120^50 km, 
whereas the distribution of larger objects is better fit by a = 3 
cumulative power law. We can give an upper size limit for the 
kink at the small-sized end around 120 km. The upper abscissa 
gives the equivalent H magnitude based on a monochromatic 
albedo distribution of pv - 0.08. 



fore smaller Plutinos is likely to be underrepresented, since the 
MPC sample suffers from an optical discovery bias. A proper 
debiasing of the size and albedo distributions is planned as part 
of future work. However, we suppose that a larger number of 
small-sized Plutinos would shift the kink to smaller diameters. 
Hence, we are only able to give an upper size limit of the loca- 
tion of the kink, which is 120 km. Furthermore, we note that our 
size distribution is based on albedo measurements derived from 
objects of the size range 150-730 km, which might be different 
from the albedo distribution of smaller Plutinos. 

4.5. Correlations 

The Plutino sample was checked for correlations of the phys- 
ical parameters, orbital properties and intrinsic optical colors, 
using a modified Spearman rank correlation analysis technique. 
The technique is based on a bootstrap method and takes into 
account the uncertainties (error bars) of the parameters (data 
points), and computes the confidence interval of the correspond- 
ing correlation coefficient. Since we are taking into account data 
error bars, the results are estimates of the most probable cor- 
relation coefficient (p) accompanied by the probability P (p- 
value) of that (p) value to occur given no relationship in the 
sample, i.e. to occur by pure chance. The lower the value of P 
the higher the significance level of the correlation. We present 
the 68 % confidence interval for the correlation value (p), i.e. 
the interval which includes 68 % of the bootstrapped values. 
This confidence interval equals the canonical 1 cr interval of a 



Table 8. Results of the modified Spearman rank correlation anal- 
ysis of the Plutino sample. is the number of sample datapoints; 
(fi) and P are the most probable correlation coefficient and the 
probability of the most probable correlation coefficient to occur 
given no relationship in the sample (cf. text), respectively. The 
uncertainty of (p) denotes its 68 % confidence level. P^- is the 
significance expressed in terms of cr. We show only the results 
of parameter pairs which meet the following criteria: |(p)| > 0.3 
and P < 0.1, or which are of special interest. Parameter pairs are 
sorted by parameter type. 



Parameters 


N 




P 




physical parameters: 

d/py 17 -0.32!°;^ 


0.217 


(1.24) 


din 


7 




0.806 


(0.25) 


Pv/r] 


7 




0.870 


(0.16) 


Hid 


17 


-0 30+"-^'' 


0.234 


(1.19) 


H/pv 


17 


-0.19!Si? 


0.468 


(0.73) 


H/rj 


7 


0.08!jJi? 


0.868 


(0.17) 


orbital and physical parameters: 

e/d 17 -0.62;0;]^ 0.008 


(2.67) 


d/q 


17 


"■"^-0.17 


0.008 


(2.67) 


d/r* 


17 


0.58^55;22 


0.015 


(2.43) 


Pv/r* 


17 


10+"" 

-0.32 


0.711 


(0.37) 


color information, orbital and physical parameters 

s/r* 13 -0.58!";*^ 0.038 (2.07) 


s/d 


13 


-0 62+"- 


0.025 


(2.24) 


si H 


13 


20+°-^'' 

-0.42 


0.523 


(0.64) 


(B-R)/d 


13 


-0 50+"' 

"•-^"-0.21 


0.082 


(1.74) 



Notes. All correlation coefficients were calculated excluding Pluto. 
Correlation analyses taking rj into account are based solely on the 7 
floating-?/ fit solutions. In order to assess the effects of discovery bias, 
we include the heliocentric distance at the time of discovery r* in our 
analysis. To improve the readability we refrain from listing all correla- 
tion values of the color indices with r* and d. Instead we only show the 
correlation values of spectral slope with r* and d, which gives similar 
results. 



gaussian distribution. As an additional aid in reading the P val- 
ues we also compute its gaussian-cr equivalent significan ce Pn-- 
The details of the technique are described in Sa ntos-Saiiz et al.l 
([accepted by A&A) . Following lEfron & Tibshirani il993 ). we 
define a weak correlation by 0.3 < |{p)| < 0.6 and a strong cor- 
relation by |(p)| > 0.6. P < 0.05 indicates reasonably strong 
evidence of correlation, P < 0.025 indicates strong evidence of 
correlation and P < 0.003 describes a clear (3 cr) correlation. 

It turns out that most parameters are uncorrected. Only very 
few parameter pairs show reasonably strong evidence of correla- 
tion. Those and a few other interesting results are listed in Table 
[8] and discussed here. Other parameter pairs lack sufficient evi- 
dence for correlation or possess weak correlation values. 



We do not detect any strong and significant trend between 
diameter d and albedo pv (cf. Figure [S]). Although we can- 
not rule out the possibility of some correlation between these 
two parameters given our sample size, it is very unlikely that 
a strong correlation between the two parameters might ex- 
ist. On the other hand, both diameter and albedo seem to be 
clearly uncorrected with the rj value. 
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Diameter (km) Heliocentric Distance at the Time of Discovery (AU) 

Fig. 8. Plots of some correlations. Upper left; we find no correlation between diameter and albedo; upper right: a clear anti- 
correlation between eccentricity and diameter is the result of a detection bias; lower left: we find smaller Plutinos to be redder than 
larger ones, which contradicts what was previously suggested by Peixinho et al.. (.2004.) based on the intrinsical brightness H of a 
different Plutino sample; lower right: we find that the color of Plutinos is correlated to the heliocentric distance at the time of their 
discovery, which hints to a color bias. 



- The absolute magnitude H is not correlated with diameter, 
albedo, or t]. This shows that our sample is not biased to- 
wards large and high (or low) albedo objects, which supports 
our approach to randomly assign albedos to H magnitudes in 
Section 14.41 This finding, however, does not exclude a dis- 
covery bias on smaller diameters than those examined here 
or a discovery bias based on geometrical aspects. 

- We find a significant anti-correlation between eccentricity 
e and diameter d (cf. Figure |8]l, showing that in our sam- 
ple highly eccentric objects tend to be smaller This anti- 
correlation is precisely the opposite of the correlation be- 
tween perihelion distance q and diameter d, and also the op- 
posite of the correlation between diameter d and the helio- 
centric distance of each object at the time of it's discovery 
r*. We consider that these relationships are merely caused 
by the fact that objects on more eccentric orbits come closer 
to the Sun, which significantly improved their detectability. 
Hence, the observed relation between an object's diameter 
and its eccentricity (and perihelion distance) is very likely to 
be caused by a discovery bias. 

- The correlation between diameter d and the heliocentric dis- 
tance of each object at the time of its discovery r*, compared 
to the lack of correlation between albedo pv and r*, shows 
that the likelihood of detection solely depends on size, and 
is rather independent of the object's albedo. This follows be- 
cause the brightness of an object scales as d^pv'- the mea- 
sured range of diameters translates into a brightness change 
factor of 24, whereas the total effect of the albedo can only 
account for a factor of 6. We can rule out a bias towards high 
(and low) albedos as a result of the nature of discovery of our 



sample targets. This supports the representativeness of the 
measured albedo distribution for the Plutino population and 
shows that the nature of the discovery bias is size-dependent. 

- Most color indices and the spectral slope are anti-correlated 
with heliocentric distance at the time of an object's discov- 
ery r* (Figure [8]l. The farther an object is, the less red it 
seems to be. This points to a probable color bias in our sam- 
ple, maybe induced by the common use of R band filters and 
the improved sensitivity in red bands of state-of-the-art de- 
tectors, which are used in TNO surveys. However, we find a 
lack of bluer objects at shorter heliocentric distances, which 
should be detectable, despite their color This might suggest 
that such objects do not exist, at least wi thin our sample. 

- We ca n n ot confirm a trend fou nd by iHainaut & Delsantil 
(2002) and 'Pei xinho et all ([2004') of bluer Plutinos being in- 
trinsically fainter using the spectral slope i as a measure of 
color The correlation analysis of B-R and H leads to similar 
results: (p) = 0.25;^q43 and P = 0.420. Our sample suggests 
precisely the opposite, with a trend between spectral slope 
s and the diameter d showing sm aller Plutinos to be re dder 
(Figure[8]l. We note however, that lPeixinho et alj (|2004|) find 
their trend due to a 'cluster' of blue Plutinos with Hr > 7.5, 
i.e. // > 8, objects at magnitude/size ranges we do not have 
in our sample. Hence, both effects might be the result of se- 
lection effects. 

It is also interesting to ask if the presence of water ice on 
our targets (see Table [3]) is correlated with albedo or diameter 
Only 6 of our targets (other than Pluto) are known to definitely 
have ice, while two definitely do not. Because we do not have a 
quantitative and consistently defined measure of how much ice 
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is present on these objects, a formal correlation analysis is not 
possible. Qualitatively speaking, we can say that 5 of the objects 
with ice have high albedos (> 0.11), with 1999 TC36 being the 
only icy object with a typical (0.07) Plutino albedo. Thus there 
is qualitative evidence that icy Plutinos have higher albedos than 
is typical. 1996 TP66 and 2000 GN171 ai-e the only objects with 
good spectra indicating the lack of water ice. Their albedos (0.07 
and 0.22, respectively) span most of the range of our measured 
albedos. This seems to indicate that Plutinos lacking water ice 
can have almost any albedo (although the lack of a detection of 
70 jum emission from 1996 TP66 by either Spitzer or Herschel 
casts considerable doubt on the accuracy of its albedo determina- 
tion). It seems remarkable that these two objects, which are free 
of ices, are also the two smallest objects (~ 150 km) of our sam- 
ple. However, the small number of spectroscopically examined 
objects does not allow for a conclusion, whether the presence of 
ices is correlated to the object diameter. 

5. Summary 

The analysis of the diameters and albedos of 18 Plutinos using 
PACS photometry, leads us to the following conclusions: 

- The diameter and albedo range of our Plutino sample yields 
150 - 730 km and 0.04 - 0.28, respectively. Excluding 
Pluto, the weighted mean of the Plutino albedo distribution is 
0.08 ± 0.03 and agrees with the average albedos observed for 
the scattered disk. Centaur and JFC population. This agree- 
ment is compatible with the idea that Plutinos are a source 
population of Centaurs and JFCs. 

- The floating-// fits yield a weighted mean rj for the Plutino 
sample of l.ll^n j n, which agrees w ith the canonical rj - 
1.20 ± 0.35 (iStansberrv et al.ll2008l) utilized in the fixed-77 
fits. 

- We refine the size estimates of the components of 
1999 TC36: cIai = 272+[^ km, cIai = 251+[^ km and 
(Jb = 132^^ km; we estimate the bulk density of the multiple 
system tohe g - 0.64;^q[j g cm"^. 

- Pluto is the clear outlier of the Plutino population, being 
by far the largest object with the brightest surface. We have 
shown, that canonical simple thermal modeling of Pluto data 
leads to inadequate results, probably mainly due to its multi- 
component nature, thermal lightcurve and atmosphere. 

- Using our measured albedos we calibrated the Plutino size 
scale for the first time and determined a cumulative size dis- 
tribution of the known Plutinos. We find that intermediate 
sized Plutinos (120 km < D < 400 km) follow a cumulative 
power-law distribution with q - 2, whereas the distribution 
of larger objects is better described by ^ = 3. We are able to 
give an upper size limit for the kink at the small-sized end of 
the distribution, which is 120 km. 

- We find no correlations between albedo and diameter, as well 
as H, and diameter, albedo and rj, respectively. This shows 
that our sample is not biased towards large and high (or low) 
albedo objects. Furthermore, we find a correlation between 
diameter and the heliocentric distance at the time of discov- 
ery, but not between the albedo and the heliocentric distance 
at the time of discovery. This leads us to the conclusion that 
the nature of the discovery bias is mainly size-dependent. 

- A significant correlation between diameter and eccentricity 
(and perihelion distance) is very likely to be caused by a de- 
tectional bias based on geometrical aspects. 

- We find hints to color biases in our sample: Plutinos, which 
have been farther from Sun at the time of their discovery 



seem to be bluer and smaller sized Plu tinos tend to be redder, 
which contradicts previous finding by'Peixinho et al.' (2004|). 
- There is qualitative evidence that icy Plutinos have higher 
albedos than the average of the sample. We are not able to 
conclude on a correlation between the diameter and the pres- 
ence of ices. 
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